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^ ; Abstract 
I ■ 

On i An aspect of the synchronization dynamics is investigated in this work. We 

as ■ 

I argue analytically and confirm numerically that the chaotic dynamics on the 

I synchronization manifold exhibits unstable dimension variability. Unstable 

o : 

, dimension variability is a cause of severe modeling difficulty for physical phe- 

H I nomena, since trajectories obtained from the mathematical model may not 

> ■ 

^ , be related to trajectories of the actual system. We present an example of 

, unstable dimension variability occurring in a system of two coupled chaotic 

maps, considering the dynamics in the synchronization manifold and its corre- 
sponding transversal direction, where a tongue-like structure is formed. The 
unstable dimension variability is revealed in the statistical distribution of the 
finite-time transversal Lyapunov exponent, having both negative and positive 
values. 
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I. INTRODUCTION 



We usually classify dynamical systems into deterministic and stochastic ones. Deter- 
ministic systems are characterized by a set of prescribed mathematical rules which evolve 
the dynamical variables in time. For these kinds of models we can numerically generate 
trajectories over long periods of time. Stochastic systems are based on some sort of random 
process. They may occur when extrinsic noise is present. A process also appears to be 
stochastic when a large number of degrees of freedom is involved. For these systems, only 
statistical information can be extracted, like averages or fluctuations of physical quantities. 

In this paper we consider a new example of a third, recently discovered, kind of dynamical 
system {pseudo- deterministic), which - in spite of being deterministic - yields only statistical 
relevant information. The problem with those systems is not related with the correctness 
and exactitude of the dynamical equations, but rather with a mathematical pathology called 
unstable dimension variability [|l],0]. 

The goodness of deterministic models is determined by some well-known paradigms: (i) 
the model must be based on sound theoretical framework, e.g., correctly applied physical 
laws; and (ii) the trajectories should reproduce correctly in some sense the actual behavior 
observed in Nature. This motivated the introduction of the model shadow ahility concept : 
let A and B be two closely related deterministic models of a physical system, but with some 
difference, which could be related to a small change in one of the system parameter values, 
or a slightly different external influence on each model, or a different noise level. The latter 
cause is restricted to arbitrarily small time dependent and bounded perturbations, which 
excludes Gaussian white noise, for example. 

We say that model shadowability occurs if trajectories of model A stay close to some 
trajectories of model B. This is necessary, but not sufficient, for either model to reproduce 
and predict correctly the time evolution of the system which the model is intended to 
describe. In other words, if there is no model shadowability neither A or B would generate 
trajectories that are physically realized, since if no trajectory of A is close to any trajectory 
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of B, it is unlikely that either model would give a trajectory that stays close to any real 
trajectory produced by Nature. 

The difficulties that obstruct model shadowability have been divided into three classes: 
minor, moderate and severe 0. Minor modeling difficulties occur for hyperbolic chaotic 
systems, since they present sensitive dependence on initial conditions. In terms of our 
example, if A and B are hyperbolic chaotic systems, trajectories of A can always be closely 
followed, or shadowed, by trajectories of B for an infinite time 

A chaotic set is hyperbolic if, at each point of a trajectory on this set, the local phase 
space can be split into unstable and stable subspaces, the angle between them being bounded 
away from zero. The unstable (stable) subspace evolves into the unstable (stable) one 
along the trajectory. We will call unstable (stable) dimension the dimension of the unstable 
(stable) subspace. However, chaotic non-hyperbolic systems are much more common in 
physical applications - they may present non-hyperbolic (homoclinic) tangencies between 
the unstable and stable subspaces. For these systems we say there is a moderate modeling 
difficulty because trajectories of A are shadowed by trajectories of B for a long, yet finite 
amount of time. However, if this shadowability time is long enough, both models are still 
useful for describing physical phenomena [Q. 

Finally, pseudo-deterministic models present severe modeling difficulties since they are 
chaotic non-hyperbolic systems presenting unstable dimension variability: the unstable and 
stable subspaces along a chaotic invariant set have no tangencies, but the dimension of the 
unstable subspace varies from point to point. For this case the shadowability time is very 
short, and no useful information could be extracted from single trajectories over a reasonable 
time span, but rather statistical information based on a probability distribution |l| . 

Unstable dimension variability has been introduced by Abraham and Smale through a 
two-dimensional continuously differentiable map [^]. A common feature associated with un- 
stable dimension variability is the oscillating behavior of a finite-time Lyapunov exponent 
between negative and positive values. This occurs because typical trajectories present arbi- 
trarily long (but finite-time) segments for which the orbit on the average is repelling in one 
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dimension and other segments for which it is repelhng in two dimensions. This behavior has 
been found in a four-dimensional invertible map describing a kicked double rotor ||^. The 
usefulness of statistical predictions has been argued through computation of the average 
energy and higher moments for this model [0]. 

Recently a non- invertible two-dimensional map was proposed as the simplest dynamical 
system exhibiting unstable dimension variability 0. Moreover it has been shown |^ that 
a lattice of diffusively coupled Henon maps present unstable dimension variability for any 
nonzero coupling strength. In this paper, we propose a new example where unstable di- 
mension variability occurs in the context of synchronization of chaotic orbits of two linear 
maps with nonlinear coupling. By changing variables to the synchronization manifold of 
the coupled system, we get a map that was first studied in the context of super-persistent 
chaotic transients and crises Pj. We have identified the mechanism which brings unstable 
dimension variability to the chaotic invariant set of this system, namely a saddle-repeller 
bifurcation which was formerly related to the boundary crisis mechanism 0. It produces a 
structure composed of super-narrow tongues through which trajectories on a chaotic saddle 
may escape after very long transients before they are reinjected, since the map is taken to be 
modulo some number. A chaotic saddle is an invariant compact set C that is both attract- 
ing and repelling and it contains a chaotic trajectory which is dense in C. This structure 



is similar to that observed in the context of the so-called riddling bifurcations |T^. The 
fluctuating behavior of the transversal finite-time Lyapunov exponent is described for this 
example and statistical information about its distribution is presented. 

This paper is organized as follows: in the second Section we introduce the coupled chaotic 
map system and analyze its synchronization manifold and the corresponding non-invertible 
map on a torus. The following two Sections are devoted to the description of the saddle- 
repeller bifurcation, the average transient lifetime and the birth of unstable-dimension vari- 
ability. Section V deals with the distribution of finite-time transversal Lyapunov exponent 
and the relative fraction of its positive values. The final Section contains our conclusions. 
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II. COUPLED CHAOTIC MAPS 



Coupled dynamical systems are susceptible to the synchronization of their trajectories, 
by which they undergo closely related motions, even when they are chaotic. In the latter 
case, even if two identical systems are started with different initial conditions, if the coupling 



is strong enough their states are asymptotically equal as the time evolves |rT|. This is a quite 
different behavior, compared to that expected from uncoupled yet identical chaotic systems, 
since if they are started with approximately equal but different initial conditions, sensitive 
dependence will cause the two systems to have completely uncorrelated motion after some 
time m. 

We consider two piecewise linear maps of the form Xn+i = kxn (mod 1), where k > 1. 
For almost all trajectories of each map the (infinite time) Lyapunov exponent A = Infc is 
positive. By almost all we mean that all orbits generated from this map are chaotic, except 



for a set of zero Lebesgue measure containing countably infinite periodic points [jT2[. We 
can write k = 1 + a/2, which is greater than one provided a > 0, and change the variable 
range from [0, 1) to [0, 7r-\/2), without altering the chaotic nature of the map orbits. In this 
case, X may be regarded as an angle variable. 

Let us introduce a nonlinear coupling of two such chaotic maps, which is equivalent to a 
nonlinear map on a torus T^: 



Un+1= {l + '^]un+U{Un,Wn) (mod VTV^), (l) 



= ( 1 + ^ ) w„ + W(m„, w„) (mod 7rv/2), (2) 



where U and W are given by 



U{Un,Wn) = ( 1 - 1^ + ^(^n - W^f + COS j , (3) 

W(m„, W„) = ( 1 - M„ - ^(Wn - Wnf - COS ^"" ^ , (4) 



where /3 > 0. 



5 



For coupled dynamical systems like this one, we can obtain the synchronization state, 
which is given by Un = Wn. Geometrically, this state defines a synchronization manifold 
S, which is a one-dimensional subset of the phase space = [0,7rv^) x [0,7rv^). If we 
represent this torus on a square with periodic boundaries, the synchronization manifold is 
a straight line with unitary slope (Fig. 1). 

In order to describe the dynamics in the synchronization manifold and in the direction 
transversal to it, we make a coordinate transformation, equivalent to a clockwise rotation 
of the axes through an angle of 7r/4: 

u + w u-w 

""^^ '^=^- 

There results a two-dimensional non-invertible map on a torus 

On+i = 26^ (mod 27r), (6) 
Zn+i = aZn + zl + (3 COS 6'„ (mod 7tV2), (7) 

where — vr/ -\/2 < z < +n/\/2. Here the synchronization manifold is simply the axis z = 0. If 
(3 = 0, the map decouples into two independent maps in 6 and z, so that an initial condition 
in the synchronization manifold will generate a chaotic orbit {^n}^o "^i^h Zn = for all 
times, that will never escape from S. Hence, S is an invariant manifold only for (3 = 0. 

An invariant manifold is typically related to the existence of some kind of symmetry in 
the system, so we may call (3 a symmetry-breaking psnametei. For nonzero (3 a chaotic orbit 
of the system is not restricted to the synchronization manifold, and can occupy a larger 
portion of the phase-space along the transversal direction z. 

The map (H-0) was introduced, in a slightly different form, P] to describe a kind of crises 
characterized by long-lived, or super-persistent, chaotic transients. They have used the z— 
part of the map without the modulo vry^ prescription, i.e., z was allowed to have any real 
value in (— oo,+oo). For some parameter values, as a = 0.7 and (3 = 0.02, the map (H-0) 
shows a chaotic attractor near the z = line. It consists of a series of curved strips with 
self-similar structure, as Henon-like bands (Fig. 2) |13|. A crude argument that can be used 



to justify the chaotic nature of this attractor consists on hnearizing the map (|^) about z = 
using z = P(, where C is a small quantity. We obtain the map Cn+i = Q^Cn + cos6'„, that 
would give a chaotic attractor for \a\ < 1, as proved by Kaplan and Yorke 

If the modulo vry^ prescription were absent from z-map, the system would present, in 
addition, a non-chaotic attractor at infinity {z = +00), and the basin of attraction of the 
two attractors would present a fractal boundary p| . In our case, however, there is no other 
attractor except the one near z = 0. Besides this chaotic attractor, there is another chaotic 
set which is a chaotic saddle (a non-attracting chaotic invariant set). It is located where the 
fractal basin boundary would be if the modulo requirement were disregarded. Trajectories 
which would accelerate to large- 2; values, are reinjected into the negative- 2; region due to the 
modulo prescription. 

III. SADDLE-REPELLER BIFURCATION 

The mechanism whereby the chaotic attractor of the map loses hyperbolicity 

through unstable dimension variability is basically an unstable-unstable pair bifurcation. 
As a result, the chaotic attractor may collide with the chaotic saddle and disappears into 
a larger chaotic saddle from which trajectories may escape, through a complex structure of 
super-narrow tongues. It is also called saddle- repeller bifurcation , and has been found as 
the cause of other strange behavior in chaotic systems, like riddling of basins of attraction 
||T0| and boundary crises 

A linear stability analysis can show the basic features of this transition. The period-1 
fixed points of the map are 

^=0, z = ^(l-a± ^(r^^^^j^^) . (8) 

Defining 

1 — a 2 /n\ 

these fixed points are written as 

7 



(e = 0,2b = z, + ^/3. and {e = {),Zc = z,-^l3,-l3). (10) 

Let us fix our attention to the case depicted in Fig. 2, i.e., a = 0.7 and /3 = 0.02 < = 
0.0225. It turns out that Zc hes in the upper point (i.e., the point with the highest value of 
z) of the chaotic attractor, whereas z^ is in the lowest point of the chaotic saddle. 
The Jacobian matrix of the map is 



f 2 ^ 



(11) 



y — /5 sin 9n a + 2Zn j 

whose eigenvalues are = 2 and ^ = « + Izn- So, the ^-direction is always unstable, as 
it should be due to the existence of the chaotic attractor. The eigenvalue related to the 
transversal direction, evaluated at the fixed points (p!0|), gives 



6(zfe,c) = l±2v//3.-/3, (12) 

so that, for /5 < z;, is a repeller, since it has an unstable dimension of dimension 2, and 
no stable subspace at all. z^ is a saddle point, since both stable and unstable subspaces have 
dimensions equal to one (Fig. 3(a)). 

If /3 = the fixed points Zc and Z}, coalesce Qk> z = z^ (Fig. 3(b)). Since the eigenvalue 
is equal to one at this point, linear stability analysis fails to determine its stability, and we 
have a saddle-repeller bifurcation with eigenvalue +1 at /5 = For (3 > jS^, the fixed points 
no longer exist, and we have a different dynamical behavior (Fig. 3(c)). 

When the fixed points coalesce at /? = the chaotic attractor collides with the chaotic 
saddle, and it becomes a chaotic transient through an interior crisis (since the attractor 
has collided with an unstable periodic orbit) 0. Trajectories arising from initial conditions 
belonging to the former basin of the chaotic attractor will typically approach its remnant, 
that is now a part of a larger chaotic saddle and will bounce around it in an irregular fashion. 
However, after some (typically very long) time this trajectory will stay near the region where 
the fixed points have coalesced, and will leave the chaotic saddle, being eventually reinjected 
to the vicinity of the saddle due to the modulo prescription. 
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At the location where the fixed points coalesced ( tongue opens-up, allowing 

trajectories near the chaotic saddle to escape for /5 > (i^. Simultaneously, each preimage of 
2;* also develops a tongue. Since these preimages are dense in the chaotic saddle, an infinite 
number of these tongues open-up simultaneously when /3 = Actually, these tongues will 
develop at those points where Q = 27im/2^, with m and k positive integers. For our map, 
however, these tongues are very narrow, since their width decrease geometrically, and are 
extremely difficult to find numerically, with exception of the main tongue opened-up at 2;*. 



IV. BEHAVIOR OF CHAOTIC TRANSIENTS 

To understand how chaotic transients are formed after the saddle-repeller bifurcation, 
let us consider the z-part of the map (|^) at = (from where the main tongue grows up) 

Zn+l = + zl + P, (13) 

whose fixed points Zb and Zc are the intersections between the parabolic function and the 
first bisector Zn+i = Zn (Fig. 4(a)). As (3 approaches its critical value these points 
approach each other and eventually coalesce when (3 = f3^, at z = z^:, where the map function 
is tangent to the 45°-line (Fig. 4(b)). For (3 > (3^ the parabola has moved upwards and 
does not intercept the bisector, leaving no fixed points. However, provided (3 is not much 
far from (3^:, a narrow channel forms between the parabola and the 45°-line, through which 
passes the trajectory resulting from the map iterations, staying there a very long time, 
and eventually escaping to large z-values (Fig. 4(c)). Due to the modulo requirement, the 
map iterates are reinjected and enter again in the channel. This is the basic mechanism of 
Type-I Intermittency but since the ^-direction is unstable, the slow motion through 



the channel does not imply a laminar behavior. Rather, it is related to a chaotic transient, 
that is characterized by irregular wandering of the trajectory over a limited 6'-range. 

The chaotic transient decays when the trajectory enters a tongue, escaping toward larger 
z-values, before it is reinjected. The main tongue is formed about the 6* = line, so let us 
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focus our attention on this spot. For (3 values greater than a tongue intercepts the 
synchronization manifold z = in an aperture of width I near 6' = 0, through which the 
trajectories can escape. Once having entered the aperture, the trajectory stays in it a 
number T of iterates before leaving the region y < Uc- Although yc must be less than 7r/-\/2, 
its exact value does not affect our results in a significant way. From now on, we set yc = 2.0. 
The trajectory, however, does not immediately stays close to ^ = since the ^-direction is 
unstable (its eigenvalue is ^2 = 2 > 1, pulling back the trajectory into the vicinity of the 
chaotic saddle). 

Let i be the distance between the orbit and = at a given time. After T iterates 
this distance increases to {^2)^ ^- However, since the 6'-excursion is bounded, we expect 
that (^2)^^ < /^i, where ki is a 0{1) constant. Now we consider what happens with a large 
number of trajectories arising from randomly chosen initial conditions. Due to the ergodicity 
of the ^-motion, we expect that almost any initial condition 9q will generate an orbit which 
has a uniform distribution over [0, 27r) for large times. Hence, the probability of O^+i falling 
into the aperture [— £/2, +£/2], if On is not in this interval, is equal to the interval width 
^. So, we can estimate the average lifetime < r > of the chaotic transient as the inverse 
of this probability. Taking our previous estimate as an upper bound for the distance along 
6'-direction at time T, we have 

< r >= K^i^f = 1^2^, (14) 

where k = 1/ki. 

In order to compute T, we assume that, for (3 slightly greater than the aperture width 
a. is very small, so that we may approximate 6' by and use Eq. ([1^) again. The difference 
Zn+i — Zn has a local minimum at z^, so we describe the dynamics within the narrow channel 
by using 6n = Zn — z^. This difference evolves with time according to the map 

5n+l=6n + 6^^ + {(3-[3,), (15) 

and since 6n+i is very close to 6n we can approximate the difference 6n+i — 5„ as a differential 
and write a differential equation for 6, now a continuous function of time r: 
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d5_ 

dr 



(16) 



which can be integrated from — 7r/-\/2 to +7r/^/2, to give an estimate for the time T it takes 



for the trajectory to escape once it has fallen in the aperture: 



T 



2 



arctan 



(_Kh/2_\ 



(17) 



Inserting Eq. (P!7| ) into Eq. (|T^), and taking logarithms we have 
logio < >= logio + 2 log^o - /9*)~"^^^ arctan ( ^ 




(/3 - /?*) 



1/2 



(18) 



To check the validity of the hypotheses made in the above derivation, we have made 
a numerical experiment, choosing Ng^ initial conditions randomly distributed over [0,27r), 
and computing the transient exit time once a given trajectory crosses the line y = Vc = 2.0. 
The average transient lifetime < r > was then computed for many values of the difference 



prediction above is the fitting parameter. 

The fitting is asymptotic though, since numerical results are best fitted by Eq. (p!8D 
when we approach i.e., in the neighborhood of the saddle-repeller bifurcation. The 
lifetime of the transients can be very large, for instance up to 10^ transients, even though 
we are as far from (3^ as 44%. This occurs because the aperture width i is so narrow that it 
becomes extremely difficult for an orbit to escape. Therefore, a better agreement between 
the theoretical prediction (|r^) and the numerical result only occurs for f3 closer to f]^,. Far 
from the theoretical estimate we made ceases to be valid, but remains useful as a lower 
bound for the exit time of transients. 



We have seen that, at the saddle-repeller bifurcation /3 = an infinite number of points, 
dense on the chaotic attractor, become repellers (unstable dimension two), and an infinite 
number of super-narrow tongues crop up as result of the collision between the attractor and 



'2 

, the results being depicted in Fig. 5, where the parameter k in the theoretical 



V. FINITE-TIME LYAPUNOV EXPONENTS 
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the chaotic saddle. However, there remains a dense set of saddle points (unstable dimension 
one) in the invariant set, and since these two different sets are densely intertwined, unstable 
dimension variability does occur in the chaotic invariant set for a large range of parameter 
values of the map. 

Yet, another signature of unstable dimension variability, as stated in the Introduction, 
is the fluctuating behavior of the finite-time Lyapunov exponents between negative and 
positive values. Consider a d-dimensional map x f(x), where x is a d-dimensional vector 
and f is a d-dimensional vector field. Let n be a positive integer and let Df"(xo) denote the 
Jacobian matrix of f" (the n-times iterated map function) evaluated at the point xq. The 
eigenvalues of the Jacobian matrix Df'^(xo) are 

cri(xo, n) > (T2(xo, n)> ...> crrf(xo, n) > 0. (19) 

We define the k-th time-n Lyapunov exponent associated with the initial condition xq as 

Afc(xo,n) = ^ln||Df"(xo)ufc||, (20) 

where is the eigenvector corresponding to the eigenvalue 0"^. Note that the usual infinite- 
time Lyapunov exponent 

\k = lim Afc(xo,n) (21) 

has the same value for almost every initial point xq with respect to the Lebesgue measure 
in the basin of attraction. 

For our map on a 2-torus, x„ = (^„, Zn)'^ , the Jacobian matrix of the n-times iterated 
map has eigenvalues ai = (^i)*^ = 2", and 

n 

^2 = 11" + 2^i, (22) 

i=l 

SO that the second (transverse) finite-time Lyapunov exponent is given by 

1 " 

X2{^o,n) = -J2\n\a + 2zi\, (23) 
^ i=i 
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the first exponent being simply In 2. 

The possibility of fluctuations between positive and negative values for this exponent 
makes useful to define a distribution function for it. Let P{X2{'XQ,n),n) denote the proba- 
bility density function of the second time-n Lyapunov exponent, when xq is chosen at random 
according to the Lebesgue measure. In other words, P(A2(xo, n), n)(iA2 is the probability 
that the exponent value lies between A2 and A2 + dX2. If F{X2) is any function of the time-n 
Lyapunov exponent, its average over the invariant measure of the attractor is given by 

(F(A2(xo, n)) = J F(A2(xo, n))P(A2(xo, n), n)dX2 (24) 

To obtain the distribution P{X2) numerically, we picked many randomly chosen initial 
conditions uniformly distributed in [0,2tt), and iterate each initial condition xq a large 
number of times. Every n = 10 steps we compute the time- 10 exponent according to Eq. 



(p3|). Actually, we use the recurrency of dynamics and follow a single trajectory a large 
number of steps, say 2 million iterates. The time n = 10 exponents are then computed from 
2 X 10^ consecutive and non-overlapping length-10 sections of the trajectory. From these 
exponents we compute a frequency histogram with convenient normalization for that 

/+00 
P(A2(xo,n),n)rfA2 = l. (25) 
-00 

In Figure 6(a), we show a distribution for a = 0.7 and (3 = 0.04, which is ~ 78% away 
from the critical value = 0.0225. In this case we can observe a distribution which looks 
like a Gaussian, but with asymmetric tails. The negative tail is sharply cut off, whereas 
the positive tail decreases smoothly. Only 0.23% of the second finite-time exponents are 
positive, indicating that almost all trajectory sections are transversally contracting. This 
is consistent with the trajectory behavior in the narrow channel that occurs near but 
the noteworthy feature here is the relatively small number of positive exponents. Figure 
6(b) depicts the same situation, but for /3 = 0.07, which is about three times the previous 
deviation away from the critical value. This time the Gaussian character of the distribution 
is apparent only for negative exponents, while there is a long, non-Gaussian fiat tail of 
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positive exponents (9% of their total number). The maximum of the distribution, however, 
appears not to have moved either toward less negative or positive values of A2(10), having 
approximately the same value of —0.35 in both cases. The fate of the probability distribution 
as we increase further the symmetry-breaking parameter (3 is illustrated in Fig. 6(c), where 
we used (3 = 0.15. The same general characteristics of the previous figure are still here, even 
though we are now very far (almost six times) from the critical value. The negative peak 
still exists at approximately the same location, but the distribution has broadened in that 
place. The novel aspect is the emergence of a second peak in the positive tail. We note that 
the fraction of positive exponents has increased to 37%. 
The fraction of positive time-n exponents 

f{n) = / P(A2(xo,n),n)rfA2 (26) 

JO 

has been computed for various values of /3, the results being depicted in Fig. 7, showing a 
monotonic increase of this fraction, indicating that for (3 ~ 0.17 about 40% of the exponents 
are positive. However, this number increases in the map (j§-0) due to the emergence of 
the second peak in the positive tail of the distribution. The shape of the curve in Fig. 7 
strongly suggests a kind of integrated probability distribution. We have thus computed the 
cumulative histogram 

P{X„n)d\, = l- P{X„n)dX„ (27) 

-OO J \2 

where we have used Eq. (^5]). In this case Q{X2 — oo) = and Q{X2 +oo) = 1. In 
Fig. 8 we show a cumulative histogram related to the distribution depicted in Fig. 6(c), 
i.e., for (3 = 0.15. In fact, if we compare it with the previous figure, qualitatively they are 
very similar, since for positive exponents it deviates from an integral of a Gaussian shaped 
function. Using Eq. (p^), it is easy to show that f{n) = 1 — Q{0, n). 
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VI. CONCLUSIONS 



In summary, this work proposes a theoretical mechanism for the existence of unstable 
dimension variability in synchronized coupled chaotic systems, which is confirmed by analysis 
and numerical computation of finite-time Lyapunov exponents. 

For invariant sets of a dynamical system, unstable dimension variability can be a strong 
obstacle to mathematical modeling of physical phenomena, since there is little probability 
that a real chaotic trajectory is shadowed for moderately long times by a trajectory from a 
model. This puts some serious doubts on the deterministic model itself, not because it is 
intrinsically bad, but rather because it presents a pathology that prevents adequate model 
shadowability. The consequence is that, although the model is deterministic, we expect to 
make only relevant statistical predictions, like averages or fluctuations, based on it. 
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FIGURE CAPTIONS 



FIGURE 1: Synchronization manifold and its transversal direction for a 2-torus. 

FIGURE 2: Phase portrait of the map for a = 0.7 and (3 = 0.02 < /3* = 0.0225. 
The dark region contains points which are driven to higher z— values but are eventually 
reinjected to negative 2;- values due to the modulo requirement. The chaotic saddle is the 
boundary between the dark and light grey region. The chaotic attractor is embedded in the 
light grey region. 

FIGURE 3: Fixed points for the map and their stability for (a) P < P^, (h) P = P^, 
and (c) P > P^. 

FIGURE 4: Return map for Eq. ® with ^ = and a = 0.7 for (a) P = 0.01 < p^ = 
0.0225, (b) = = 0.0225, and (c) p = 0.04 > p^. 

FIGURE 5: Base- 10 logarithm of the average exit time of transients < r > of the map 
(^-|^) with a = 0.7 and = 0.0225 as a function of {P — P^)'^^"^ . Circles: numerical 
experiment, full line: theoretical prediction given by Eq. (|T8|) with a fitting parameter 
fi; = 0.0045. 

FIGURE 6: Probability distribution for the transverse time-10 Lyapunov exponents for 
a = 0.7 and (a) P = 0.04; (b) P = 0.07; and (c) p = 0.15. 

FIGURE 7: Fraction of positive transverse time-10 Lyapunov exponents f{n = 10) as a 
function of P, for a = 0.7 and /5* = 0.0225. 

FIGURE 8: Cumulative histogram for the transverse time-10 Lyapunov exponents for 
a = 0.7 and p = 0.15. 
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